LHT1/MAC7 contributes to proper alternative splicing under long-term heat stress and mediates variation in the heat tolerance of Arabidopsis.

Natural genetic variation has facilitated the identification of genes underlying complex traits such as stress tolerances. We here evaluated the long-term (L-) heat tolerance (37°C for 5 days) of 174 Arabidopsis thaliana accessions and short-term (S-) heat tolerance (42°C, 50 min) of 88 accessions and found extensive variation, respectively. Interestingly, L-heat-tolerant accessions are not necessarily S-heat tolerant, suggesting that the tolerance mechanisms are different. To elucidate the mechanisms underlying the variation, we performed a chromosomal mapping using the F2 progeny of a cross between Ms-0 (a hypersensitive accession) and Col-0 (a tolerant accession) and found a single locus responsible for the difference in L-heat tolerance between them, which we named Long-term Heat Tolerance 1 (LHT1). LHT1 is identical to MAC7, which encodes a putative RNA helicase involved in mRNA splicing as a component of the MOS4 complex. We found one amino acid deletion in LHT1 of Ms-0 that causes a loss of function. Arabidopsis mutants of other core components of the MOS4 complex-mos4-2, cdc5-1, mac3a mac3b, and prl1 prl2-also showed hypersensitivity to L-heat stress, suggesting that the MOS4 complex plays an important role in L-heat stress responses. L-heat stress induced mRNA processing-related genes and compromised alternative splicing. Loss of LHT1 function caused genome-wide detrimental splicing events, which are thought to produce nonfunctional mRNAs that include retained introns under L-heat stress. These findings suggest that maintaining proper alternative splicing under L-heat stress is important in the heat tolerance of A. thaliana.

Expression profiles in knock-down transgenic plants of high and low diversified duplicate genes in Arabidopsis thaliana.

Duplicated genes show various degrees of functional diversification in plants. We previously identified 1,052 pairs of high diversified duplicates (HDDs) and 600 pairs of low diversified duplicates (LDDs) in Arabidopsis thaliana. Single knock-down of HDDs induced abnormal phenotypic changes because the other gene copy could not compensate for the knock-down effect, while single knock-down of LDDs did not induce abnormal phenotypic changes because of functional compensation by the copy gene. Here, focusing on one pair each of HDDs and LDDs, we performed transcriptome analyses in single-knock-down transgenic plants. The numbers of differentially expressed genes in single-knock-down transgenic plants were not different between HDDs and LDDs. Thus, functional compensation inferred by transcriptomics was similar between HDDs and LDDs. However, the trend of differentially expressed genes was similar in the pair of LDDs, while expression profiles were dissimilar in the pair of HDDs. This result indicates that a pair of LDDs tends to share similar functions but a pair of HDDs tends to have undergone functional divergence. Taking these findings together, as the reason for no phenotypic changes in single knock-down of LDDs but phenotypic changes in double knock-down of LDDs, we concluded that phenotypic changes of LDDs were induced by decreasing gene dosage.

Unearth of small open reading frames (sORFs) in drought stress transcriptome of Oryza sativa subsp. indica.

Drought is a major abiotic stress that influences rice production. Although the transcriptomic data of rice against drought is widely available, the regulation of small open reading frames (sORFs) in response to drought stress in rice is yet to be investigated. Different levels of drought stress have different regulatory mechanisms in plants. In this study, drought stress was imposed on four-leaf stage rice, divided into two treatments, 40% and 30% soil moisture content (SMC). The RNAs of the samples were extracted, followed by the RNA sequencing analysis on their sORF expression changes under 40%_SMC and 30%_SMC, and lastly, the expression was validated through NanoString. A total of 122 and 143 sORFs were differentially expressed (DE) in 40%_SMC and 30%_SMC, respectively. In 40%_SMC, 69 sORFs out of 696 (9%) DEGs were found to be upregulated. On the other hand, 69 sORFs out of 449 DEGs (11%) were significantly downregulated. The trend seemed to be higher in 30%_SMC, where 112 (12%) sORFs were found to be upregulated from 928 significantly upregulated DEGs. However, only 8% (31 sORFs out of 385 DEGs) sORFs were downregulated in 30%_SMC. Among the identified sORFs, 110 sORFs with high similarity to rice proteome in the PsORF database were detected in 40%_SMC, while 126 were detected in 30%_SMC. The Gene Ontology (GO) enrichment analysis of DE sORFs revealed their involvement in defense-related biological processes, such as defense response, response to biotic stimulus, and cellular homeostasis, whereas enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways indicated that DE sORFs were associated with tryptophan and phenylalanine metabolisms. Several DE sORFs were identified, including the top five sORFs (OsisORF_3394, OsisORF_0050, OsisORF_3007, OsisORF_6407, and OsisORF_7805), which have yet to be characterised. Since these sORFs were responsive to drought stress, they might hold significant potential as targets for future climate-resilient rice development.

A de novo gene originating from the mitochondria controls floral transition in Arabidopsis thaliana.

De novo genes created in the plant mitochondrial genome have frequently been transferred into the nuclear genome via intergenomic gene transfer events. Therefore, plant mitochondria might be a source of de novo genes in the nuclear genome. However, the functions of de novo genes originating from mitochondria and the evolutionary fate remain unclear. Here, we revealed that an Arabidopsis thaliana specific small coding gene derived from the mitochondrial genome regulates floral transition. We previously identified 49 candidate de novo genes that induce abnormal morphological changes on overexpression. We focused on a candidate gene derived from the mitochondrial genome (sORF2146) that encodes 66 amino acids. Comparative genomic analyses indicated that the mitochondrial sORF2146 emerged in the Brassica lineage as a de novo gene. The nuclear sORF2146 emerged following an intergenomic gene transfer event in the A. thaliana after the divergence between Arabidopsis and Capsella. Although the nuclear and mitochondrial sORF2146 sequences are the same in A. thaliana, only the nuclear sORF2146 is transcribed. The nuclear sORF2146 product is localized in mitochondria, which may be associated with the pseudogenization of the mitochondrial sORF2146. To functionally characterize the nuclear sORF2146, we performed a transcriptomic analysis of transgenic plants overexpressing the nuclear sORF2146. Flowering transition-related genes were highly regulated in the transgenic plants. Subsequent phenotypic analyses demonstrated that the overexpression and knockdown of sORF2146 in transgenic plants resulted in delayed and early flowering, respectively. These findings suggest that a lineage-specific de novo gene derived from mitochondria has an important regulatory effect on floral transition.

Enhanced growth rate under elevated CO2 conditions was observed for transgenic lines of genes identified by intraspecific variation analyses in Arabidopsis thaliana.

KEY MESSAGE: Using the whole genome and growth data of Arabidopsis thaliana ecotypes, we identified two genes associated with enhancement of the growth rate in response to elevated CO2 conditions. Improving plant growth under elevated CO2 conditions may contribute to enhanced agricultural yield under future global climate change. In this study, we examined the genes implicated in the enhancement of growth rates under elevated CO2 conditions by analyzing the growth rates of Arabidopsis thaliana ecotypes originating from various latitudes and altitudes throughout the world. We also performed a genome-wide association study and a transcriptome study to identify single nucleic polymorphisms that were correlated with the relative growth rate (RGR) under elevated CO2 conditions or with CO2 response of RGR. We then selected 43 candidate genes and generated their overexpression and/or RNA interference (RNAi) transgenic mutants for screening. After screening, we have found that RNAi lines of AT3G4000 and AT5G50900 showed significantly higher growth rates under the elevated CO2 condition. As per our findings, we conclude that natural variation includes genetic variation associated with the enhancement of plant productivity under elevated CO2 conditions.

Positive selective sweeps of epigenetic mutations regulating specialized metabolites in plants.

DNA methylation is an important factor regulating gene expression in organisms. However, whether DNA methylation plays a key role in adaptive evolution is unknown. Here, we show evidence of naturally selected DNA methylation in Arabidopsis thaliana In comparison with single nucleotide polymorphisms, three types of methylation-methylated CGs (mCGs), mCHGs, and mCHHs-contributed highly to variable gene expression levels among an A thaliana population. Such variably expressed genes largely affect a large variation of specialized metabolic quantities. Among the three types of methylations, only mCGs located in promoter regions of genes associated with specialized metabolites show a selective sweep signature in the A thaliana population. Thus, naturally selected mCGs appear to be key mutations that cause the expressional diversity associated with specialized metabolites during plant evolution.

Genome-Wide Association Study of Local Thai Indica Rice Seedlings Exposed to Excessive Iron.

Excess soluble iron in acidic soil is an unfavorable environment that can reduce rice production. To better understand the tolerance mechanism and identify genetic loci associated with iron toxicity (FT) tolerance in a highly diverse indica Thai rice population, a genome-wide association study (GWAS) was performed using genotyping by sequencing and six phenotypic data (leaf bronzing score (LBS), chlorophyll content, shoot height, root length, shoot biomass, and root dry weight) under both normal and FT conditions. LBS showed a high negative correlation with the ratio of chlorophyll content and shoot biomass, indicating the FT-tolerant accessions can regulate cellular homeostasis when encountering stress. Sixteen significant single nucleotide polymorphisms (SNPs) were identified by association mapping. Validation of candidate SNP using other FT-tolerant accessions revealed that SNP:2_21262165 might be associated with tolerance to FT; therefore, it could be used for SNP marker development. Among the candidate genes controlling FT tolerance, RAR1 encodes an innate immune responsive protein that links to cellular redox homeostasis via interacting with abiotic stress-responsive Hsp90. Future research may apply the knowledge obtained from this study in the molecular breeding program to develop FT-tolerant rice varieties.

Degree of Functional Divergence in Duplicates Is Associated with Distinct Roles in Plant Evolution.

Gene duplication is a major mechanism to create new genes. After gene duplication, some duplicated genes undergo functionalization, whereas others largely maintain redundant functions. Duplicated genes comprise various degrees of functional diversification in plants. However, the evolutionary fate of high and low diversified duplicates is unclear at genomic scale. To infer high and low diversified duplicates in Arabidopsis thaliana genome, we generated a prediction method for predicting whether a pair of duplicate genes was subjected to high or low diversification based on the phenotypes of knock-out mutants. Among 4,017 pairs of recently duplicated A. thaliana genes, 1,052 and 600 are high and low diversified duplicate pairs, respectively. The predictions were validated based on the phenotypes of generated knock-down transgenic plants. We determined that the high diversified duplicates resulting from tandem duplications tend to have lineage-specific functions, whereas the low diversified duplicates produced by whole-genome duplications are related to essential signaling pathways. To assess the evolutionary impact of high and low diversified duplicates in closely related species, we compared the retention rates and selection pressures on the orthologs of A. thaliana duplicates in two closely related species. Interestingly, high diversified duplicates resulting from tandem duplications tend to be retained in multiple lineages under positive selection. Low diversified duplicates by whole-genome duplications tend to be retained in multiple lineages under purifying selection. Taken together, the functional diversities determined by different duplication mechanisms had distinct effects on plant evolution.

Effect of small coding genes on the circadian rhythms under elevated CO2 conditions in plants.

Increase in atmospheric carbon dioxide (CO2) has a significant effect on plant growth and development. To explore the elevated-CO2 response, we generated transcriptional profiles over a time course (2 h-14 days) of exposure to elevated CO2 in Arabidopsis thaliana. Genes related to photosynthesis were down-regulated and circadian rhythm-related genes were abnormally regulated in the early to middle phase of elevated CO2 exposure. To understand the novel mechanism of elevated CO2 signaling, we focused on 42 unknown small coding genes that showed differential expression patterns under elevated CO2 conditions. Four transgenic plants overexpressing the small coding gene exhibited a growth-defective phenotype under elevated CO2 but not under current CO2. Transcriptome analysis showed that circadian rhythm-related genes were commonly regulated in four transgenic plants. These circadian rhythm-related genes were transcribed in the dark when CO2 concentrations in the leaf was high. Taken together, our identified four small coding genes are likely to participate in elevated CO2 signaling to the circadian rhythm.

Identification of endogenous small peptides involved in rice immunity through transcriptomics- and proteomics-based screening.

Small signalling peptides, generated from larger protein precursors, are important components to orchestrate various plant processes such as development and immune responses. However, small signalling peptides involved in plant immunity remain largely unknown. Here, we developed a pipeline using transcriptomics- and proteomics-based screening to identify putative precursors of small signalling peptides: small secreted proteins (SSPs) in rice, induced by rice blast fungus Magnaporthe oryzae and its elicitor, chitin. We identified 236 SSPs including members of two known small signalling peptide families, namely rapid alkalinization factors and phytosulfokines, as well as many other protein families that are known to be involved in immunity, such as proteinase inhibitors and pathogenesis-related protein families. We also isolated 52 unannotated SSPs and among them, we found one gene which we named immune response peptide (IRP) that appeared to encode the precursor of a small signalling peptide regulating rice immunity. In rice suspension cells, the expression of IRP was induced by bacterial peptidoglycan and fungal chitin. Overexpression of IRP enhanced the expression of a defence gene, PAL1 and induced the activation of the MAPKs in rice suspension cells. Moreover, the IRP protein level increased in suspension cell medium after chitin treatment. Collectively, we established a simple and efficient pipeline to discover SSP candidates that probably play important roles in rice immunity and identified 52 unannotated SSPs that may be useful for further elucidation of rice immunity. Our method can be applied to identify SSPs that are involved not only in immunity but also in other plant functions.

Contribution of Functional Divergence Through Copy Number Variations to the Inter-Species and Intra-Species Diversity in Specialized Metabolites.

There is considerable diversity in the specialized metabolites within a single plant species (intra-species) and among different plant species (inter-species). The functional divergence associated with gene duplications largely contributes to the inter-species diversity in the specialized metabolites, whereas the intra-species diversity is due to gene dosage changes via gene duplications [i.e., copy number variants (CNVs)] at the intra-species level of evolution. This is because CNVs are thought to undergo associated with less functional divergence at the intra-species level of evolution. However, functional divergence caused by CNVs may induce specialized metabolite diversity at the intra-species and inter-species levels of evolution. We herein discuss the functional divergence of CNVs in metabolic quantitative trait genes (mQTGs). We focused on 5,654 previously identified mQTGs in 270 Arabidopsis thaliana accessions. The ratio of nonsynonymous to synonymous variations tends to be higher for mQTGs with CNVs than for mQTGs without CNVs within A. thaliana accessions, suggesting that CNVs are responsible for the functional divergence among mQTGs at the intra-species level of evolution. To evaluate the contribution of CNVs to inter-species diversity, we calculated the ratio of nonsynonymous to synonymous substitutions in the Arabidopsis lineage. The ratio tends to be higher for the mQTGs with CNVs than for the mQTGs without CNVs. Additionally, we determined that mQTGs with CNVs are subject to positive selection in the Arabidopsis lineage. Our data suggest that CNVs are closely related to functional divergence contributing to adaptations via the production of diverse specialized metabolites at the intra-species and inter-species levels of evolution.

Tyr82 Amino Acid Mutation in PB1 Polymerase Induces an Influenza Virus Mutator Phenotype.

In various positive-sense single-stranded RNA viruses, a low-fidelity viral RNA-dependent RNA polymerase (RdRp) confers attenuated phenotypes by increasing the mutation frequency. We report a negative-sense single-stranded RNA virus RdRp mutant strain with a mutator phenotype. Based on structural data of RdRp, rational targeting of key residues, and screening of fidelity variants, we isolated a novel low-fidelity mutator strain of influenza virus that harbors a Tyr82-to-Cys (Y82C) single-amino-acid substitution in the PB1 polymerase subunit. The purified PB1-Y82C polymerase indeed showed an increased frequency of misincorporation compared with the wild-type PB1 in an in vitro biochemical assay. To further investigate the effects of position 82 on PB1 polymerase fidelity, we substituted various amino acids at this position. As a result, we isolated various novel mutators other than PB1-Y82C with higher mutation frequencies. The structural model of influenza virus polymerase complex suggested that the Tyr82 residue, which is located at the nucleoside triphosphate entrance tunnel, may influence a fidelity checkpoint. Interestingly, although the PB1-Y82C variant replicated with wild-type PB1-like kinetics in tissue culture, the 50% lethal dose of the PB1-Y82C mutant was 10 times lower than that of wild-type PB1 in embryonated chicken eggs. In conclusion, our data indicate that the Tyr82 residue of PB1 has a crucial role in regulating polymerase fidelity of influenza virus and is closely related to attenuated pathogenic phenotypes in vivoIMPORTANCE Influenza A virus rapidly acquires antigenic changes and antiviral drug resistance, which limit the effectiveness of vaccines and drug treatments, primarily owing to its high rate of evolution. Virus populations formed by quasispecies can contain resistance mutations even before a selective pressure is applied. To study the effects of the viral mutation spectrum and quasispecies, high- and low-fidelity variants have been isolated for several RNA viruses. Here, we report the discovery of a low-fidelity RdRp variant of influenza A virus that contains a substitution at Tyr82 in PB1. Viruses containing the PB1-Y82C substitution showed growth kinetics and viral RNA synthesis levels similar to those of the wild-type virus in cell culture; however, they had significantly attenuated phenotypes in a chicken egg infection experiment. These data demonstrated that decreased RdRp fidelity attenuates influenza A virus in vivo, which is a desirable feature for the development of safer live attenuated vaccine candidates.

Hormone-like peptides and small coding genes in plant stress signaling and development.

Recent works have shed light on the long-distance interorgan signaling by which hormone-like peptides precisely regulate physiological effects in a manner similar to phytohormones. Many such peptides have already been identified in the primary model plant, Arabidopsis thaliana. In addition, Arabidopsis genome reanalysis revealed over 7000 novel candidate small coding genes, some of which are likely to be associated with hormone-like peptides. Hormone-like peptides have also been reported to play critical roles in interorgan communications during morphogenesis and stress responses. In this review, we focus on the functional roles of hormone-like peptides and small coding genes in cell-to-cell and/or long-distance communications during plant stress signaling and development and discuss the evolutionary conservation of these peptides among plants.

Dephosphorylation of clustered phosphoserine residues in human Grb14 by protein phosphatase 1 and its effect on insulin receptor complex formation.

The physical interaction of the human growth factor receptor-bound protein 14 (hGrb14) and the insulin receptor (IR) represses insulin signaling. With respect to the recruiting mechanism of hGrb14 to IR respond to insulin stimulus, our previous reports have suggested that phosphorylation of Ser358 , Ser362 , and Ser366 in hGrb14 by glycogen synthase kinase-3 repressed hGrb14-IR complex formation. In this study, we investigated phosphatase-mediated dephosphorylation of the hGrb14 phosphoserine residues. An in vitro phosphatase assay with hGrb14-derived synthetic phosphopeptides suggested that protein phosphatase 1 (PP1) is involved in the dephosphorylation of Ser358 and Ser362 . Furthermore, coimmunoprecipitation experiments suggested that insulin-induced hGrb14-IR complex formation was repressed by the substitution of Ser358 or Ser362 with glutamic acid. These findings suggested that phosphate groups on Ser358 and Ser362 in hGrb14 are dephosphorylated by PP1, and the dephosphorylation facilitates hGrb14-IR complex formation.

Tuning water-use efficiency and drought tolerance in wheat using abscisic acid receptors.

Water availability is a key determinant of terrestrial plant productivity. Many climate models predict that water stress will increasingly challenge agricultural yields and exacerbate projected food deficits. To ensure food security and increase agricultural efficiency, crop water productivity must be increased. Research over past decades has established that the phytohormone abscisic acid (ABA) is a central regulator of water use and directly regulates stomatal opening and transpiration. In this study, we investigated whether the water productivity of wheat could be improved by increasing its ABA sensitivity. We show that overexpression of a wheat ABA receptor increases wheat ABA sensitivity, which significantly lowers a plant's lifetime water consumption. Physiological analyses demonstrated that this water-saving trait is a consequence of reduced transpiration and a concomitant increase in photosynthetic activity, which together boost grain production per litre of water and protect productivity during water deficit. Our findings provide a general strategy for increasing water productivity that should be applicable to other crops because of the high conservation of the ABA signalling pathway.

Distinct gene expression signatures induced by viral transactivators of different HTLV-1 subgroups that confer a different risk of HAM/TSP.

BACKGROUND: Among human T cell leukemia virus type 1 (HTLV-1)-infected individuals, there is an association between HTLV-1 tax subgroups (subgroup-A or subgroup-B) and the risk of HAM/TSP in the Japanese population. To investigate the role of HTLV-1 subgroups in viral pathogenesis, we studied the functional difference in the subgroup-specific viral transcriptional regulators Tax and HBZ using microarray analysis, reporter gene assays, and evaluation of viral-host protein-protein interaction. RESULTS: (1) Transcriptional changes in Jurkat Tet-On human T-cells that express each subgroup of Tax or HBZ protein under the control of an inducible promoter revealed different target gene profiles; (2) the number of differentially regulated genes induced by HBZ was 2-3 times higher than that induced by Tax; (3) Tax and HBZ induced the expression of different classes of non-coding RNAs (ncRNAs); (4) the chemokine CXCL10, which has been proposed as a prognostic biomarker for HAM/TSP, was more efficiently induced by subgroup-A Tax (Tax-A) than subgroup-B Tax (Tax-B), in vitro as well as in unmanipulated (ex vivo) PBMCs obtained from HAM/TSP patients; (5) reporter gene assays indicated that although transient Tax expression in an HTLV-1-negative human T-cell line activated the CXCL10 gene promoter through the NF-κB pathway, there was no difference in the ability of each subgroup of Tax to activate the CXCL10 promoter; however, (6) chromatin immunoprecipitation assays showed that the ternary complex containing Tax-A is more efficiently recruited onto the promoter region of CXCL10, which contains two NF-κB binding sites, than that containing Tax-B. CONCLUSIONS: Our results indicate that different HTLV-1 subgroups are characterized by different patterns of host gene expression. Differential expression of pathogenesis-related genes by subgroup-specific Tax or HBZ may be associated with the onset of HAM/TSP.

AtPep3 is a hormone-like peptide that plays a role in the salinity stress tolerance of plants.

Peptides encoded by small coding genes play an important role in plant development, acting in a similar manner as phytohormones. Few hormone-like peptides, however, have been shown to play a role in abiotic stress tolerance. In the current study, 17 Arabidopsis genes coding for small peptides were found to be up-regulated in response to salinity stress. To identify peptides leading salinity stress tolerance, we generated transgenic Arabidopsis plants overexpressing these small coding genes and assessed survivability and root growth under salinity stress conditions. Results indicated that 4 of the 17 overexpressed genes increased salinity stress tolerance. Further studies focused on AtPROPEP3, which was the most highly up-regulated gene under salinity stress. Treatment of plants with synthetic peptides encoded by AtPROPEP3 revealed that a C-terminal peptide fragment (AtPep3) inhibited the salt-induced bleaching of chlorophyll in seedlings. Conversely, knockdown AtPROPEP3 transgenic plants exhibited a hypersensitive phenotype under salinity stress, which was complemented by the AtPep3 peptide. This functional AtPep3 peptide region overlaps with an AtPep3 elicitor peptide that is related to the immune response of plants. Functional analyses with a receptor mutant of AtPep3 revealed that AtPep3 was recognized by the PEPR1 receptor and that it functions to increase salinity stress tolerance in plants. Collectively, these data indicate that AtPep3 plays a significant role in both salinity stress tolerance and immune response in Arabidopsis.

Efficient anchoring of alien chromosome segments introgressed into bread wheat by new Leymus racemosus genome-based markers.

BACKGROUND: The tertiary gene pool of bread wheat, to which Leymus racemosus belongs, has remained underutilized due to the current limited genomic resources of the species that constitute it. Continuous enrichment of public databases with useful information regarding these species is, therefore, needed to provide insights on their genome structures and aid successful utilization of their genes to develop improved wheat cultivars for effective management of environmental stresses. RESULTS: We generated de novo DNA and mRNA sequence information of L. racemosus and developed 110 polymorphic PCR-based markers from the data, and to complement the PCR markers, DArT-seq genotyping was applied to develop additional 9990 SNP markers. Approximately 52% of all the markers enabled us to clearly genotype 22 wheat-L. racemosus chromosome introgression lines, and L. racemosus chromosome-specific markers were highly efficient in detailed characterization of the translocation and recombination lines analyzed. A further analysis revealed remarkable transferability of the PCR markers to three other important Triticeae perennial species: L. mollis, Psathyrostachys huashanica and Elymus ciliaris, indicating their suitability for characterizing wheat-alien chromosome introgressions carrying chromosomes of these genomes. CONCLUSION: The efficiency of the markers in characterizing wheat-L. racemosus chromosome introgression lines proves their reliability, and their high transferability further broadens their scope of application. This is the first report on sequencing and development of markers from L. racemosus genome and the application of DArT-seq to develop markers from a perennial wild relative of wheat, marking a paradigm shift from the seeming concentration of the technology on cultivated species. Integration of these markers with appropriate cytogenetic methods would accelerate development and characterization of wheat-alien chromosome introgression lines.

Functional divergence of duplicate genes several million years after gene duplication in Arabidopsis.

Lineage-specific duplicated genes likely contribute to the phenotypic divergence in closely related species. However, neither the frequency of duplication events nor the degree of selection pressures immediately after gene duplication is clear in the speciation process. Here, using Illumina DNA-sequencing reads from Arabidopsis halleri, which has multiple closely related species with high-quality genome assemblies (A. thaliana and A. lyrata), we succeeded in generating orthologous gene groups in Brassicaceae. The duplication frequency of retained genes in the Arabidopsis lineage was ∼10 times higher than the duplication frequency inferred by comparative genomics of Arabidopsis, poplar, rice and moss (Physcomitrella patens). The difference of duplication frequencies can be explained by a rapid decay of anciently duplicated genes. To examine the degree of selection pressure on genes duplicated in either the A. halleri-lyrata or the A. halleri lineage, we examined positive and purifying selection in the A. halleri-lyrata and A. halleri lineages throughout the ratios of nonsynonymous to synonymous substitution rates (KA/KS). Duplicate genes tended to have a higher proportion of positive selection compared with non-duplicated genes. Interestingly, we found that functional divergence of duplicated genes was accelerated several million years after gene duplication compared with immediately after gene duplication.